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1. Introduction 

We consider Euler’s equations describing the motion of a perfect incompressible fluid in vacuum: 

(1.1) {dt + V'"dk)vj =-djp, j = l,...,n in P, 

(1.2) divU = 9fcU^ = 0 in V 

where di = djdx^ and P = U o<t<r {t} x Vt, T>t C M”. Here = Vk and we use the summation 

convention over repeated upper and lower indices. The velocity vector field of the fluid is V, p is the 
pressure and is the domain the fluid occupies at time t. We also require boundary conditions on the 
free boundary dV = Uo<t<r {t} x dVt] 

(1.3) P = 0, on dV, 

(1.4) {dt + V’^dk)\oveT{dV), 

Condition (1.3) says that the pressure p vanishes outside the domain and condition (1.4) says that the 
boundary moves with the velocity V of the fluid particles at the boundary. 

Given a domain Dq C M"", that is homeomorphic to the unit ball, and initial data vq, satisfying 
the constraint (1.2), we want to find a set P = Uo<t<r {t} x Dt, C M” and a vector field v solving 

(1.1) -(1.4) with initial conditions 

(1.5) {x; (0,x) G V} = Vq, and v = vq, on {0} x Vq 

Let M be the exterior unit normal to the free surface d'Dt. Christodoulou[C2] conjectured that the 
initial value problem (1.1)-(1.5), is well posed in Sobolev spaces if 

(1.6) Vj\fp<—co<0, on dV, where Vj\f = M''d^i. 

Condition (1.6) is a natural physical condition since the pressure p has to be positive in the interior 
of the fluid. It is essential for the well posedness in Sobolev spaces. A condition related to Rayleigh- 
Taylor instability in [BHL,W1] turns out to be equivalent to (1.6), see [W2]. Taking the divergence of 

(1.1) gives: 

(1.7) —/\p={djV^)dkV\ in Vt, P = 0, on dVt 
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In the irrotational case, when cxiilvij = diVj — djVi = 0, then Ap < 0 so p > 0 and (1.6) holds by the 
strong maximum principle. Wu [W1,W2] proved well posedness locally in time, (globally in space), in 
Sobolev spaces in the irrotational case. Ebin [El] showed that the equations are ill posed when (1.6) 
is not satisfied and the pressure is negative and Ebin [E2] announced an existence result when one 
adds surface tension to the boundary condition. With Christodoulou [CL] we proved a priori bounds in 
Sobolev spaces in the general case of non vanishing curl, assuming (1.6). Usually if one has a priori 
estimates, existence follows from similar estimates for some regularization or iteration scheme for the 
equation. However, the sharp estimates in [CL] use all the symmetries of the equations and so only hold 
for perturbations of the equations that preserve the symmetries. Here we show existence in Sobolev 
spaces for the linearized equations using a new type of estimates. 

The incompressible perfect fluid is to be thought of as an idealization of a liquid. Eor small bodies 
like water drops surface tension should help holding it together and for larger denser bodies like stars 
its own gravity should play a role. Here we neglect the influence of such forces. Instead it is the 
incompressibility condition that prevents the body from expanding and it is the fact that the pressure 
is positive that prevents the body from breaking up in the interior. Let us also point out that, from 
a physical point of view one can alternatively think of the pressure as being a small positive constant 
on the boundary instead of vanishing. What makes this problem difficult is that the regularity of the 
boundary enters to highest order. Roughly speaking, the velocity tells the boundary where to move 
and the boundary is the zero set of the pressure that determines the acceleration. 

Some existence results in Sobolev spaces are known in the irrotational case, for the closely related 
water wave problem which describes the motion of the surface of the ocean under the influence of earth’s 
gravity. In that problem, the gravitational field can be considered as uniform, however this problem 
reduces to our problem by going to an accelerated frame. The domain Dt is unbounded for the water 
wave problem coinciding with a half-space in the case of still water. Nalimov[Na] and Yosihara[Y] proved 
local existence in Sobolev spaces in two space dimensions for initial conditions sufficiently close to still 
water. Beale, Hou and Lowengrab[BHL] have given an argument to show that problem is linearly well 
posed in a weak sense in Sobolev spaces, assuming a condition, which can be shown to be equivalent 
to (1.6). The condition (1.6) prevents the Rayleigh-Taylor instability from occurring when the water 
wave turns over. Finally Wu[Wl,2] proved local existence in general in two and three dimensions for 
the water wave problem. The method of proofs in these papers uses that the velocity is irrotational 
and divergence free and hence harmonic to reduce the equations to equations on the boundary only. 

The main result here is existence for the linearized equations in the case of non vanishing curl. 
The irrotational case was proved by Yosihara [Y]. The proof in [Y], see also [W1,W2], reduces the 
equation to the boundary and it does not generalize. Instead, we project the linearized equation onto 
an equation in the interior using the orthogonal projection onto divergence free vector fields in the 
inner product. This removes a difficult term, the differential of the linearization of the pressure, 
and reduces a higher order term, the linearization of the moving boundary, to a symmetric unbounded 
operator on divergence free vector fields. The linearized equation becomes an evolution equation in 
the interior for this operator, which we call the normal operator. It is basically the differential of the 
harmonic extension to the interior of the normal component. In the irrotational case it becomes the 
normal derivative which is elliptic on harmonic functions and our equation reduces to an equation on 
the boundary similar to those in [Y,W1,W2]. 

The normal operator is positive due to (1.6) and this will lead to energy bounds. However, existence 
of regular solutions does not follow from standard energy methods or semi-group methods since the 
operator is time dependent and non-elliptic in the case of non vanishing curl. Usually one gets equations 
and estimates for higher derivatives by commuting differential operators through the equation, but 
we can only use operators whose commutator with the normal operator is controlled by the normal 
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operator. Geometric arguments lead us to use Lie derivatives with respect to divergence free vector 
fields tangential at the boundary. The commutators of these with the normal operator are controlled 
by the normal operator and they preserve the divergence free condition. The same considerations apply 
to time differentiation so one should use the Lie derivative with respect to the material derivative (1.4) 
which reduces to the time derivative of the vector field in the Lagrangian coordinates. To get estimates 
for all derivatives we use the fact that we have a better evolution equation for the curl and that any 
derivative can be controlled by tangential derivatives, the curl and the divergence. 

As pointed out above, existence does not follow directly from estimates but one must have existence 
and uniform estimates for some regularizing sequence. We replace the normal operator by a sequence 
of bounded operators converging to it which are still symmetric, positive and they uniformly satisfy the 
same commutator estimates with the differential operators above. Due to the geometric construction of 
the differential operators there is a natural regularization which corresponds to replacing the boundary 
by an inhomogeneous term supported in a small neighborhood of it. 

Existence for the linearized equations or some modification will be part of any existence proof for 
the nonlinear problem. The estimates here require more regularity of the solution we linearize around 
than we get for the linearization. However, we use the techniques presented here in a forthcoming paper 
[L3], to prove existence for the nonlinear problem with the Nash-Moser technique. 

In order to formulate the linearized equations one has to introduce some parametrization of the 
boundary. Let us therefore first express Euler’s equations in the Lagrangian coordinates in which the 
boundary becomes fixed. Given a domain Dq in R"", that is diffeomorphic to the unit ball D, we can by 
a theorem in [DM] find a diffeomorphism /o : D ^ Dq that up to a constant factor is volume preserving, 
i.e. after an additional scaling det{dfo/dy) = l. Assume that V and v G C{'D) are given satisfying (1.4). 
The Lagrangian coordinates x = x{t,y) = ft{y) are given by solving 


(1-8) — = V{t,x{t,y)), x{0,y) = fo{y), yen. 

Then /t : D ^ is a volume preserving diffeomorphism, since divH = 0, and the boundary becomes 
fixed in the new y coordinates. Let us introduce the notation 


(1.9) 


Dt 




+ 1 /^ 

x=const 


d 

dx^ 




y=const 


and 


d (9y“ d 

dx^ dx^ dy°- ’ 


for the material derivative and partial differential operators expressed in the Lagrangian coordinates. 
In these coordinates Euler’s equations (1.1), the incompressibility condition (1.2) and the boundary 
condition (1.3) become 


( 1 . 10 ) 


DjX* = —dip, det (dx/dy) = 1, 


in [0, T] X D 


and 



= 0 . 


where p = p(t,y), di now is to be thought of as the differential operator in (1.9) in y and Dt is the 
time derivative. We then define V = Dtx. Note that the second equation in (1.10) follows since 
Dt In (det (dx/dy)) = divE = 0. Taking the divergence of the first equation in (1.10) gives (1.7) so p is 
determined as functional of (re, V). The initial conditions (1.5) become 


t^o = /o, Dt x = Eo 


and 


( 1 . 11 ) 

subject to the constraints, 

( 1 . 12 ) 


det {dfo/dy) = 1, 
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divEo = 0 



and Christodoulou’s physical condition become 

(1.13) ^ -Co < 0 

where M is the exterior unit normal to dT^t parametrized by x{t,y). 

Let us now derive the linearized equations of (1.10). We assume that {x{t,y),p{t,y)) is a given 
smooth solution of (1.10) satisfying (1.13) for 0 < t < T. Let x{t,y,r) and p{t,y,r) be smooth 
functions also of a parameter r, such that ix,p)\^_Q = {x,p) and set {5x,Sp) = {dx/dr,dp/dr)\^_^. 
Then the linearized equations is the requirement on {Sx,5p), that {x,p) satishes the equations (1.10) 
up to terms bounded by as r ^ 0. In other words, if 


(1.14) ^i{x,p) = DfX^ + dip, i = l,...,n, 


d>o(x,p) = det (dx/dy) - 1, 4>n+i(x,p) = p\g^, 


then then linearized operator is defined by 
(1.15) ^'{x,p){6x,6p) = 

Euler’s equations (1.10) become ^{x,p) = 


, where x = x + r5x, p = p + r6p 

I r=0 

0 and the linearized equations are 


(1.16) 


^'{x,p){6x, 5p) = 0 


Applying the operator 5f = df to (1.10), using that by (2.8) [5,di] = — {di5x^)dk-, gives the 

linearized equations 


(1.17) 


D^6x'^ — {dkp)di5x’^ = —di6p, div(5x = 0, and 



= 0 . 


where we used that (51n (det (clx/(9y)) = div(5a; see (2.6). Here 5p is determined as a functional of 
{5x,Dt5x) since taking the divergence of (1.17) gives an elliptic equation for 5p similar to (1.7). We 
now want to sow existence for (1.17) with initial data 


(1.18) 


= 5fo, Dt = (5Eo, 


satisfying the constraints 


(1.19) div(5/o = 0, div(5Vo = (9i(5/o )9fcEJ 

We remark, that the difference between (1.10) and (1.17) is the term dkpdiSx^ in (1.17). This term 
is higher order but because of the sign condition (1.6) it will contribute with a positive term to the 
energy. We also remark that the equation (1.17) also shows up in estimating energies of higher order 
derivatives for (1.10) in [CL]. In fact, the material derivative Dt corresponds to the variation 5 given 
by time translation. Our main result is: 

Theorem 1.1. Let H he the unit ball in R” and suppose that {x,p) is a smooth solution of (1-10) 
satisfying (1-13) for 0 < t < T. Suppose that {Sfo,SVo) are smooth satisfying the constraints (1.19). 
Then the linearized equations (1.17) have a smooth solution {5x,5p) for 0<t<T satisfying the initial 
conditions (1.18). LetAf be the exterior unit normal to dT>t parametrized byx{t,y) and letdxjy =J\f-5x 
be the normal component. Set 


Er{t) = \\Dt5x{t, •)||h’-(o) + ||<^x(t, OIU-CO) + 
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( 1 . 20 ) 



where H^{Q) and are the Sobolev spaces in Q respeetively on dQ. Then there are eonstants Cr 

depending only on {x,p), r and T such that 

(1.21) Er{t) < CrEr(0), for 0 <t <T, r > 0. 

Furthermore, let be the completion of divergence free veetor fields in the norm 

• Then if the constraints in (1-19) hold and 

( 1 . 22 ) i6fo,SVo)eN'^{n)xH-{n) 

it follows that (1.17)-(1.18) has a solution 

(1.23) {5x,Dt5x) e C{[0,T],N^{n) x 

As we have argued, any smooth solution of (1.1)-(1.5) with T>o diffeomorphic to the unit ball can be 
reduced to a smooth solution of (1.10) where fl is the unit ball. The term ||5x_/v'||_fff(o) is equivalent to 
the variation of the second fundamental form 6 = dAf of the free boundary dVt measured in 
so our energy is essentially 2 ( 9 ^) + This is to be compared with the a priori bounds 

for the nonlinear problem in [CL] for ||0||_f/r—A slightly more general theorem holds, 
see section 2. Let us now outline the main ideas in the proof. We will rewrite the linearized equations 
(1.17) in a geometrically invariant way and use this to obtain energy bounds and a regularization of 
the equation which will give existence. 

We have defined our functions and vector fields to be functions of the Lagrangian coordinates 
{t, y) G [0, T] X n but we can alternatively think of them as functions of the Eulerian coordinates 
(t, x) G T), and we will make this identification without explicitly saying that we compose with the 
inverse of the change of coordinate y x{t,y). The time derivative has a simple expression in the 
Lagrangian coordinates but the space derivatives have a simpler expression in the Eulerian coordinates, 
see (1.9). For the most part we will think of our functions and vector fields in the Lagrangian frame 
but we use the inner product coming from the Eulerian frame, i.e. in the Lagrangian frame we use the 
pull-back metric of the Euclidean inner product: 

(1.24) X-Z = 5i,X^Z^ = gasX'^Z\ where = = — — . 

Here A* refers to the components of the vector X in the Eulerian frame, A“ refers to the components 
in the Lagrangian frame, gab is the metric in the Lagrangian frame and 5ij is the Euclidean metric in 
the Eulerian frame. The letters a, 6 , c, d, e, /, g will refer to indices in the Lagrangian frame whereas 
the indices i, j,k,l,m,n will refer to the Eulerian frame. The norms and most of the operators we 
consider have an invariant interpretation so it does not matter in which frame they are expressed. In 
the introduction we use express the vector fields in the Eulerian frame but later we express the vector 
fields in the Lagrangian frame. The Lf inner product of vector fields is given by 

(1.25) {X,Z)= [ X-Zdx = [ X-Zdy 

Jvt Jn 

where the equality follows from the incompressibility condition: det (dx/dy) = 1 . 

We now want to derive energy bounds for the linearized equations (1.17). Let us first point out that 
the boundary condition p\q^ = 0 implies that the energy is conserved for a solution of Euler’s equations 
(1.10). We have 

(1.26) ^ [\V\‘^dx = [ Dt\Vfdx =-2 [ V^dipdx = 2 [ dwVpdx -2 [ V^pdS = 0, 

dt Jvt Jvt Jvt Jvt JdVt 
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where VV = J^iV^ is the normal component of V. In fact, the first equality follows from the incom¬ 
pressibility condition after expressing the integrals as integrals over as in (1.25), the second is Euler’s 
equations (1.10), the third follows from the divergence theorem and the last is the boundary condition 
and the divergence free condition. 

We will now use the orthogonal projection onto divergence free vector fields to rewrite the linearized 
equations (1.17) in an invariant way that can be used to derive energy bounds and for which there is a 
natural regularization. The orthogonal projection onto divergence free vector fields in the inner product 
(1.25) is given by 

(1.27) PX^ = X* — 5^^dj q, where Aq = divX, q =0. 

dDt 

We now want to project the first equation in (1.17) onto divergence free vector fields. This removes the 
right hand side di5p, since we project along gradients of functions that vanish on the boundary. The 
second term in the first equation in (1.17) can be written as —di[{dkp)6x^) + {didkp)dx^, where the last 
part is lower order and the projection of the first part turns out to be a positive symmetric operator on 
divergence free vector fields. We define the normal operator A to be 

(1.28) AX^ = P{- 5"^dj{X^dkp)) = -5^^dj{X^dkP - q). 

where q is chosen so that the divergence of AX vanishes and q vanishes on the boundary. Then yl is a 
positive symmetric operator on divergence free vector fields, if condition (1.6) holds. In fact, if X and 
Z are divergence free then 

(1.29) (X, AZ) = - [ X^di{Z^dkp) dx= [ X^Zj^ {-Vj^p) dS, X^ = NiX^ 

JT>t JdDt 

There is one more issue we have to deal with before writing up the linearized equations (1.17) in a 
more pleasant form. The time derivative Dt does not preserve the divergence free condition so we have 
to modify it so it does. The operator 

(1.30) Co,X‘ = D,X' - (4V)X‘ = 

preserves the divergence free condition if V is divergence free. This is because it is the space time Lie 
derivative with respect to the divergence free vector field Dt = (1, V) restricted to the space components. 
Another way to look at it is that it is just the time derivative of the vector field X expressed in the 
Lagrangian frame. The divergence is invariant under coordinate changes and the volume form is time 
independent so it commutes with time differentiation in the Lagrangian coordinates. 

We now project the linearized equations (1.17) and get an evolution equation on divergence free 
vector fields for the normal operator A: 

(1.31) X* + AX* = -2F((dkV^)X^) where X = Sx, X = CdJx, X = C%^5x 

Introducing the orthogonal projection onto divergence free vector fields solved to problems. First it 
turned the higher order term, the second term in (1.17) into a positive symmetric operator. Secondly it 
got rid of the third term in (1.17) which caused considerable difficulties in [CL]. In fact, the projection 
of a gradient of a function that vanishes on the boundary vanishes. The right hand side of (1.31) is 
lower order since the projection is a bounded operator. Associated with (1.31) is the energy 


(1.32) 


E(t) = (X,X) + (X, (A + /)X) 
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and one can show an energy estimate \E'{t)\ < CE(t) which gives an energy bound. We remark that 
for divergence free vector fields (1-32) is equivalent to (1.20) with r = 0. In order to show this energy 
bound we must calculate the commutator of the time derivative and the normal operator, which follows 
from the argument below. 

In order to prove the energy bound and similar energy bounds for higher derivatives one has to 
control the commutator of differential operators with the normal operator. This is however a delicate 
matter since these commutators have to be controlled by the normal operator itself and only certain 
geometric operators satisfy this. Let T be a divergence free vector held that is tangential at the boundary 
and let 

(1.33) - X’^dkT^ 

be the Lie derivative with respect to T applied to a vector held X. Then CtX is divergence free if 
X is divergence free. It turns out that the commutators between Ct and the normal operator can be 
controlled by the normal operator; 

(1.34) [£t, = {CTE^)djkAX’^ + ArpX^ 

where for / vanishing on the boundary we dehned 

(1.35) AfX = -P{E^d,{X'^dkf)). 

(1.34) follows from (1.28) using that the Lie derivative commutes with exterior differentiation and that 
the tangential derivatives Tp and Tq also vanish on the boundary since p and q do. In view of the 
physical condition (1.13) it follows from (1.29) that 

(1.36) \{X,AfX)\<C{X,AX), where C = \\V^ f/V^fp\\^9^) ■ 

Applying Ct to the linearized equations (1.31) therefore gives a similar equation for CtSx for which 
we also get energy bounds if T is a divergence free vector held that is tangential at the boundary. The 
second term in the commutator (1.34) can be controlled using (1.36). In order to control the hrst term 
in the commutator one has to use that AX can be controlled in terms of C^ 6x through the equation 
(1.31). Therefore we also have to differentiate the equation with respect to time and include time 
derivatives up to highest order in the energies. We dehne energies 

(1.37) Ej{t)= ^ {CD,CipX,CD,Ci^X) + {Ci^X,ACipX), X = 5x 

\I\<r,Ier 

where T is a family of divergence free vector helds that are tangential at the boundary and span the 
tangent space of the boundary including the time derivative Dt and Cl^ is any product of r = |/| Lie 
derivatives with respect to these. Then one can prove energy estimates EJ (t) < CEJ (0). 

The energies (1.37) only contain tangential derivatives. In order to control normal derivatives also 
we use: 

(1.38) |9Z| < (7(1 divZ| + I curlZj + E i«2i) 

ses 

where 5 is a family of vector fields that span the tangent space of the boundary and curl_Z^^. = — 

cljZj, where = 6ijZ^ is the one form corresponding to the vector field Z. The divergence of X = 
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Cd^5x vanishes and there is a better evolution equation for curlX. In fact the curl of the higher order 
operator A in (1.31) considered as an operator with values in the one forms vanishes since it is a gradient. 
For a solution of Euler’s equations (1.10) the curl is preserved: 

(1.39) jCot curlu = 0, 

where Cot is the space time Lie derivative with respect to Dt = (l,E) of the two form a: 

(1.40) = Dt (Tij + {diV )aij + {djV )aii = > 

restricted to the space components, i.e. it is the time derivative of the two form expressed in the 
Lagrangian frame. For the linearized equations we have the following identity; 

(1.41) Cot curl dz = 0, 5zi = — curlu^j = 6x^ 

Since the Lie derivative commutes with exterior differentiation cmlCfpSz is also conserved. 

The above argument gives energy bounds, assuming existence. However, existence does not follow 
directly from estimates. To show existence we must approximate the linearized equations with some 
equation for which we know there is existence and prove that we have uniform bounds for the norms 
as the approximation gets better so that we can construct a sequence that tends to a solution of the 
linearized equations. For / > 0 in and = 0 we define the smoothed out normal operator by 

(1.42) A^X^ = P{ - Xe{d)P^d,{fd-^X'^dkd)) = P{x'Ad)S^^idjd)fd-^X^dkd) 

where d = d{y) = dist {y,dQ,) and Xe{d) = x{d/£)- Here x is a smooth cut off function, x(s) = 1, when 
•s > 1) x('S) = 0, when s < 0 and x^(s) > 0. Then is a positive symmetric operator on divergence 
free vector fields, if condition (1.6) holds. In fact, if X and Z are divergence free then 

(1.43) {X,A%Z) = - [ Xe{d)X^di{fd-^Z’^dkd)dx = [{X^did){Z'^dkd) Xeid)fd-Ux. 

Jvt Jvt 

It follows that is symmetric and positive and satisfies the same commutator properties as Af and 
the curl of vanishes when d > £. Furthermore is a bounded operator, i.e. ||HjX||^ < Cer-HXH^. 

We will actually first obtain energy estimates for the linearized equations with vanishing initial data 
and an inhomogeneous divergence free term that vanishes to any order as t ^ 0: 

(1.44) X* + AX^ + 2P{{dkV^)X^) = 5$, = 0, k < r, div(54> = 0, 

of the form 

(1.45) Ej{t)<Cr f \\6^J dr, where ||5$||^ = ^ 

\I\<r,IeT 

One can reduce to this situation by subtracting a power series solution in time to (1.31). (1.44) with 
A replaced by A^ = A^ is just an ordinary differential equation in P[^{kl) so existence for this equation 
follows. Because A^ uniformly satisfies the same commutator estimates as A we will obtain uniform 
energy bounds and will be able to pass to the limit as e ^ 0 and obtain a solution for (1.44). The reason 
we have to hrst subtract off the initial conditions in this way is that the energy (1.37) contains time 
derivatives up to highest order and these would have to be obtained from the equation. The operator 
A^ is smoothing but only in the tangential directions and in the normal directions it is worse than A so 
if we had replaced A by A^ directly in (1.31) the higher order initial conditions would have depended on 
in an uncontrollable way. As described above, we will first prove the energy bounds in such a way 
that we can obtain the same uniform bounds for the smoothed out equation and pass to the limit as 
e ^ 0 to obtain existence. Once we have existence we can then obtain the more natural energy bounds 
for the initial value problem in Theorem 1.1. 



2. Lagrangian coordinates, the linearized equation and statement of the theorem. 

Let US introduce Lagrangian coordinates in which the boundary becomes fixed. Let be a domain 
in R” and let /o : ^ Do be a difFeomorphism that is volume preserving; det(9/o/9|/) = 1. For 

simplicity we will assume that Vol(Do) is the volume of the unit ball in R^. By a theorem of [DM] 
we can prescribe the volume form up to a constant for any mapping of one domain into another so we 
may assume that D is the unit ball. Assume that v{t,x) and p{t,x), {t,x) £ V are given satisfying the 
boundary conditions (1.3)-(1.4). The Lagrangian coordinates x = x{t,y) = ftiy) are given by solving 


( 2 . 1 ) 


dx/dt = V{t,x{t,y)), x{0,y) = fo{y), y£^ 


Then /t : D ^ Dt is a volume preserving diffeomorphism, since divD = 0, and the boundary becomes 
fixed in the new y coordinates. Let us introduce the notation 


( 2 . 2 ) 


"‘ = 1 


y=constant 


A 

dt 


jrk 9 

+ - . 

x=constant dx^ 


for the material derivative and 


(2.3) 


d dy°' d 

' 


In these coordinates Euler’s equation (1.1), the incompressibility condition (1.2) and the boundary 
condition (1.3) become 


(2.4) Djx* = -dip, K = det (dx/dy) = 1, p\g^ = 0 

where x = x(t,y), p=p(t,y). The initial conditions (1.5) become 

(2.5) a;|j^o = /o, Dtx\^^^ = Vo. 

In fact, recall that Df det (M) = det (M) tr DtM), for any matrix M depending on t so 

(2.6) Dt det {dx/dy) = det (dx/dy) (dy^/dx'^){dDtX^/dy°‘) = diDtX^ = divD^x = divE = 0. 


Note that p is uniquely determined as a functional of x by (2.4)-(2.5). In fact taking the divergence of 
Euler’s equations (2.4) using (2.6) gives Ap = —{diDtx^){djDtx'^). 

Let 5 be a variation with respect to some parameter r, in the Lagrangian coordinates; 


(2.7) 


S = d/dr 


We think of x{t, y, r) and p{t, y, r) as depending on r and differentiate with respect to r. Differentiating 
(2.3) using the formula for the derivative of the inverse of a matrix, 5M~^ = —M~^{SM)M~^, gives 


( 2 . 8 ) 


[5, Si] = -{di6x'")dk. 


Differentiating (2.4), using (2.8) and (2.6) with Df replaced by 6 gives the linearized equations; 


D/5x^ — {dkp)di5x^ = —di6p, 


(2.9) 
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div(5x = 0, ^P\dn ~ 



It is however better to use the fact that v and p are solutions of Euler’s equations, DtVi = —dip, to 
arrive at the following equation 


( 2 . 10 ) - di[{dkp)6x^) = -di5p + {dkDtVi)5x^ 


divdx = 0, ^P\dn ~ 


We will now transform the vector field 6x to Lagrangian coordinates, because in these coordinates 
the time derivative preserves the divergence free condition. Let 


( 2 . 11 ) 


1E“ = 5x‘ 


dyO. 

(9X* ' 


dx* = W‘ 


,clx* 

dy^ 


q = 6p. 


The letters a, b, c, d, e, / will refer to quantities in the Lagrangian frame whereas the letters i, j, k,l,m,n 
will refer to ones in Eulerian frame, e.g. da = 9/<9y“ and di = djdx^. With this convention we have 


( 2 . 12 ) 


d. = ^d d = —d- 

'-'i o ,• '-'a-: '-'a n'-'I- 

ax* ay“ 


Multiplying the first equation in (2.10) by dx^ jdy^" and summing over i gives 
(2.13) 


6ij—D^5x^ - da{{dcp)W^) = daq + 


since {dkp)Sx^ = {dcp)W‘^ and {dkDtVi)5x^ = {dcDtVi)W^. On the other hand 


(2.14) 

(2.15) 




and 






dyb Qyb .. Qyb 

Multiplying (2.15) by dx'^ jdy^, summing over i, and substituting into (2.13) gives 


(2.16) 

where 

(2.17) 




9ab — d 


(9x* dx^ 

a Qya Qyb 


is the metric 5ij expressed in the Lagrangian coordinates. Let be the inverse of the metric gab, 

dx^ dx^ dx'^ dx^ 

(2.18) gab = Dtgab = ^ ~ 

be the time derivative of the metric and the vorticity in the Lagrangian coordinates. Expression (2.16) 
becomes 


(2.19) 


gabD^.W'’ - da{{daP)W^) = -daq - {gac ” 0Jac)DtW‘^. 
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(2.19) can alternatively be expressed, using the inverse of gab-, in the form 

( 2 . 20 ) - g^'>db{{daP)W<^) = -g^^^dtq - - cobc)DtW^. 

The divergence is invariant under coordinate changes so the second condition in (2.10) is 

(2.21) divlT = K~^da{K.W°') = 0, where k = det {dx/dy) = 1 
Finally, the last equation in (2.10) is, since q = 5p, 

( 2 . 22 ) = 0 

Then linearized equations are now the requirement that (2.20), (2.21) and (2.22) hold and we want to 
find (IT, q) satisfying these equations and the initial conditions 

(2.23) = Wo, = Wi, where divlTo = divlTi = 0, IT = DfW 

We can however express (2.20)-(2.23) in as one equation as follows. First we note that q = Sp is 
determined as a functional of IT and DtW. In fact, it follows from (2.21) that divU^ IT = 0 so taking 
the divergence of (2.20) using (2.22) gives us an elliptic equation for q: 

(2.24) Ag = da{ng'^''dbq) = da{Kg‘"''db{{dcp)W‘^) - Kg‘^’’{gbc - Wbc) , q\g^ = 0. 

We now write q = qi + q 2 + qs, where |= 0 and Aqi is equal to each of the three terms in the right 
hand side of (2.24). The equations (2.20)-(2.22) can then be written as one equation, LiIT = 0 where 

(2.25) LiW = W + AW+ GW-CW 
and 

(2.26) AW‘^ = -g‘^^db{{dap)W^-qr), Agi = A((9eP)IT^) • Qi\dn = ^ 

(2.27) GIT“ = 5“'(ff(,cIT^ + g2), Ags = q2\g^ = 0 

(2.28) CW‘^ = g‘^\iObcW^-q3), Aqs = daig'^^cvbcW^^) <?3|af, = 0 

We will prove the following theorem: 

Theorem 2.1. Suppose that x,p € C'“([0,T] x fl), p\qq = 0, < —cq < 0 and divDtx = 0. 

Suppose that F G C'°“([0,T] x Q) and ITq, ITi G C'°“(fl) are all divergenee free. Then 

(2.29) L,W = F, wl^^ = Wo, IT|^^q = ITi, 

where Li be given by (2.25)-(2.28), has a divergence free solution IT G C'“([0, T] x 17). 

Let be the Sobolev spaces and let A’’(I7) be the completion ofC°°{Ll) divergenee free vector 

fields in the norm ||IT||j;i-r-(-f 7 ) + WW^WH^iaci), where ITw = IT • A is the normal component. Then if 

(2.30) (ITo,ITi) G iV’'(I7) X 77’'(f7), F G ([0,T],i7’'(f7)) 
are all divergence free it follows that (2.29) have a a divergence free solution 

(2.31) (IT, IT) G C([0,r],iV’'(f7) x i7’'(I7)). 

Moreover, with a constant C depending only on the norm of x and p and the constant cq we have 

(2.32) \\W{t)\\H^ + ||IT(t)||;v^ < C(||IT(0)||^^. + ||IT(0)||w^ + \\F{t)\\h^ dr). 

Remark. The restrictions that divT = 0 and divF = 0 can be removed and in order to use the Nash- 
Moser technique one indeed needs to show that the linearized operator is invertible away from a solution 
and outside the divergence free class. In [L3] the techniques presented here are used to show this. 
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3. The projection onto divergence free vector fields and the normal operator. 
Let P be the orthogonal projection onto divergence free vector fields in the inner product 


{W,U) = / 
Jn 


(3.1) 

Then the projection P 

(3.2) PU^^ = U<^-g^^^dbq, Ag = (?) = divt/= 5,(1/“), q =0. 

dfl 

That this is the orthogonal projection follows since gab9^‘^ = 

(3.3) {W,{I-P)U) = - [ gabW^g^^d,qdy= [ {daW^)qdy-[ NaW^qdS = 0, if 5,1T“ = 0 

J fi J J 

where A, is the exterior unit conormal and dS is the surface measure. The projection of a gradient of 
a function that vanishes on the boundary vanishes; 


(3.4) 


P{9^'’d,q) = 0 , 


if 


^\dn 


The projection has norm one: 

(3.5) \\PU\\ < ||[/||, 11(7 - P)U\\ < ||[/||, IIITII = (IT, 1T)V2 

The projection is continuous on the Sobolev spaces 77’’(17) if the metric is sufficiently regular; 

(3.6) \\PU\\H^in)<Cr\\U\\H^^n), 
since it is just a matter of solving the Dirichlet problem: 

(3-7) \\q\\H-+^n) < Cr\\U\\H-(n), r > 0, if Aq = divU, = 0. 

For r > 1 this is the standard estimate for the Dirichlet problem. For r = 0 this is obtained by 
multiplying by q, using that the right hand side is in divergence form, integrating by parts and using 
that = 0. Furthermore if the metric also depends smoothly on time t then 

k k 

(3.8) ^||77^P[/||H.(f,) <a,fc j;i|77^'t/|k^(D). 

j=o i=o 

This follows by induction in k from commuting through time derivatives in (3.2); 

m-l , 

(3.9) ADTq = -Yl f" 5a((77rVVi.77^'9)+5,(77rt/“), = « 

j=o / 7 

which using (3.7) gives \\D^q\\H-+i(ji) < Cr,mYJj=o \\PT9\\H-+^(n) + Am EJIg II 
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For functions / vanishing on the boundary we define operators on divergence free vector helds 


(3.10) = P{-g‘^^db{{dJ)W-)), 

Af is symmetric, i.e. {U,AfW) = {AfU, IF), since for U and IF divergence free it follows from (3.3) 

(3.11) {U,AfW) = - [ U^da{{dJ)W‘^)Kdy= [ {-VNfWNWNndS, Un = NaU^ 

Jn Jon 

If p is the pressure in Euler’s equations then normal operator A in (2.26) is 


(3.12) 


A = Ap>{), i.e. (IF, ^IF) > 0, if V/vP 


dn 


< 0 


which is true by our assumption (1.6). It follows from Cauchy Schwartz inequality that 

(3.13) \{U,AfpW)\ < {U,Ay\pU)^/\W,A\j\pWf/^ < ||/|U»(ao)(t/,^C/)'/"(lF,^IF)'/^ 

since V/v(E) = fV^p on the boundary. The positivity properties (3.12) and (3.13) are of fundamental 
importance to us. In particular, since p vanishes on the boundary so does p = Dtp and therefore 


(3.14) 


A = Ap satisfies \{W, Aw)\ < \\VNp/VNp\\L^idn){W, AW) 


A is the time derivative of the operator A, considered as an operator with values in the one forms. 

It follows from (3.10) and (3.5) that ||^/1F|| < ||9^/||L“(r2)||11^|| + ||5 /||l“(o)|| 511^||. However, Af 
acting on divergence free vector fields by (3.11) depends only on V/v/|gf^, i.e. A^ = Af if = 

Vv/laf^. We can therefore replace / by the Taylor expansion of order one in the distance to the 
boundary in polar coordinates multiplied by a smooth function that is one close to the boundary and 
vanishes close to the origin. It follows that 

(3.15) \\AfW\\<C'£\\^NSf\\L^^on)\\W\\+C\\VNf\\L^idn){\\dW\\ + \\W\\). 

ses 

where 5 is a set of vector helds that span the tangent space of the boundary, see section 6. 

For two forms a we dehne bounded projected multiplication operators given by 


(3.16) 


M„IF“ = P(5“VIF^), I|M„IF|| < ||a||L»(f7) 


In particular the operators in (2.27) and (2.28) are bounded projected multiplication operators; 


(3.17) 


G = Mg, C = M^, G=Mg 


where g is the metric, lv the vorticity and g the time derivative of the metric. 
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4. The lowest order energy estimate. 

Since det {dx/dy) = 1 it follows from introducing Lagrangian coordinates, that for a function / 

(4.1) f fdx= f fdy, so ^ f f dx = f Dtfdx 

JVt Jn Jvt Jvt 

We note that if i; is a solution if Euler’s equations, DfVi = —diP, and p vanish on the boundary then 

(4.2) ^ f \V\‘^dx = 2 f V^DtVidx = —2 f V^dipdx = 2 f {divV)pdx — 2 f VNpdS = 0 

dt Jvt Jvt Jvt JdDt 

We now want to obtain energy estimates for the linearized equations 

(4.3) LiW = W + AW+ GW-CW = F 

where A, G and G are as in section 3 and F is divergence free. Because of the unbounded but positive 
and symmetric operator A there is an additional term in the energy: 

(4.4) E = E{W) = {W,W) + {W,{A + I)W) 
where the inner product is given by (3.1). 

Since (W, W) = j^QabW'^WUy and Dt{gabW<^W^) = QabW'^W'^ + 2g,^W<^DtW\ we have 

(4.5) ^{W, W) = 2{W, DtW) + (IT, GW) 

where G is given by (3.17). By (3.3) and (3.11) {W,AW) = — Jq W°'da((dcp)W^) dy, and 

(4.6) Dt{W‘^da{{dcp)W‘^))) = W^da{{dcp)W‘^)) + W^daiidcPW^^) + W‘^da{{d,Dtp)W‘^). 

Since A is symmetric we get 

(4.7) ^{W,AW) = 2{W,AW) + {W,AW) 
where AW^ = ApW^ is given by (3.10) with f = p = Dtp. Hence 

(4.8) j^E{W) = 2{W, W + AW + W) + (IT, GIT) + (IT, ilT) + (IT, GIT) 

= 2(1T, LilT) + 2(1T, IT) - (IT, GVT) + (VT, AW) + (VT, GIT). 

where we used that (IT, GIT) vanishes since G is antisymmetric. The operator G is bounded by (3.16)- 
(3.17) and |(1T, yilT)| is bounded by (3.14) so 

(4.9) \E\ < (^1+ + \\VNDtp/VNp\\L'=°{du)jE+ 2'/E\\F\\. 

With n(t) = 1 + ||ff||L~(o) + \\'^NDtp/VNP\\L^{dn) and Eq = \/E we hence have 

(4.10) Eo{t) <ed>‘^^(^Eo{0) + 
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5. Turning the initial conditions into an inhomogeneous divergence free term. 


As explained in the introduction we want to reduce the initial value problem 

(5.1) LiW = W+ AW+ GW-CW = F, = Wo, = Wi 

to the case of vanishing initial conditions and an inhomogeneous term F that vanishes to any order as 
t ^ 0. This is achieved by subtracting off a power series solution in t to (5.1): 

F “ 1 “ 2 

( 5 . 2 ) WSr{t,y) = Y.^-^W:{y) 

s=0 

We note that if Ws are divergence free it follows that ITor is divergence free. Here ITo and ITi are the 
initial conditions, IT 2 is obtained form the equation (5.1) at t = 0: IT 2 = F — AWq — GWi + GW\. 
Similarly, one gets higher order terms by first differentiating the equation with respect to time. It is 
clear that doing so we obtain an expression = Mj.{W ,+ D^F and from this we 

inductively define Wk +2 = Mk{Wo ,..., Wk+i)\^^Q + Here is some linear operator of order 

at most one and that is all we need to know. However, we are going to calculate the explicit form of 
Mk since we will do similar calculations later on for other operators and this is a simple model case. 

Now it turns out that its easier to differentiate the corresponding operator with values in one forms; 

(5.3) L^Wa = QabLlW^ = gabW^ - a4(9eP)fH") + daq + {gab - COabW^ = 9 abF^ 

where q is chosen so the last terms are divergence free, and afterwards project the result to the divergence 
free vector fields. Let 


(5.4) q^=Dtq = D^p, gi, = Dlg,,, = DtcOab, F, = D^F 

In general it follows from applying to (5.3), restricting to t = 0 gives that 


(5-5) XI 0 (sirwl+j - o,({a,p’-)wi)) + a.f + ^ Q( 9 ;; 

s=0 


s=0 


s=0 


— LO 


ab 


We now want to project each term onto divergence free vector fields. Let 

(5.6) A^W^^ = P{- g^^^iidaP^W ^)), GsW^ = P{g^'^gl^W '^), G^W^^ = P(g‘^^oj:,W^) 

We obtain 


r— 1 r 

(5.7) W ,+2 = - (;)G,_.W ,+2 - (D {Gr-s+lWs + l - Gr-sWs + l + - Gr-sFs) 

s=0 s=0 

This inductively defines Wr+2 from Wq, ..., HV+i. With ITor given by (5.2) we have hence achieved that 

(5.8) Dt{LiWor-F)l^o = 0, for s < r, Wor\^^^ = Wo, Wor\^^^ = W^ 

Replacing IT by IT — ITor and F by F — LilTor- hence reduces (5.1) to the case of vanishing initial data 
and an inhomogeneous term that vanishes to any order r as t ^ 0. 

We also note that if the initial data are smooth then we can construct a smooth approximate 
solution IT that satishes the equation to all orders as t ^ 0. This is obtained by multiplying the 
term in (5.2) by a smooth cutoff xit/^k), to be chosen below, and summing up the infinite series. Here 
y is smooth y(s) = 1 for |s| < 1/2 and y(s) = 0 for |s| > 1. The sequence > 0 can then be chosen 
small enough so that the series converges in C'™'([0, T], F"*) for any m if we take (||lTfc||fc + 1)^^ < 1/2. 
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6. Construction of the tangential vector fields. 


Let us now construct the tangential divergence free vector fields, that are time independent expressed 
in the Lagrangian coordinates, i.e. that commute with Dt'. 

(6.1) [Dt,T]=0. 

This means that in the Lagrangian coordinates they are of the form T°'{y)d/dy^ and since det (dx/dy) = 
1 the divergence free condition is just 

( 6 . 2 ) daT^ = 0 . 

Since is the unit ball in R” the vector fields can be explicitly given. The vector fields 

(6.3) y^djdy^-y^djdy- 

corresponding to rotations, span the tangent space of the boundary and are divergence free in the 
interior. Furthermore they span the tangent space of the level sets of the distance function from the 
boundary in the Lagrangian coordinates 

(6.4) d{y) =dist{y,dn) = l-\y\ 

away from the origin y / 0. We will denote this set of vector fields by Sq We also construct a set 
of divergence free vector fields that span the full tangent space at distance d{y) > do and that are 
compactly supported in the interior at a fixed distance do/2 from the boundary. The basic one is 

(6.5) h{y ^,..., y*") (^fiy^)g'iy^)d/dy^ - f'iy^)giy^)d/dy ^'^, 

which is divergence free. Furthermore we can choose /, g, h such that it is equal to d/dy^ when |y*| < 1/4, 
for i = 1, ...,n and so that it is 0 when |y*| > 1/2 for some i. In fact let / and g be smooth functions 
such that f{s) = 1 when |s| < 1/4 and f{s) = 0 when |s| > 1/2 and g'{s) = 1 when |s| < 1/4 and 
g{s) = 0 when |s| > 1/2. Finally let h{y^, ...,y^) = f{y^) • • • By scaling, translation and rotation 

of these vector fields we can obviously construct a finite set of vector fields that span the tangent space 
when d> do and are compactly supported in the set where d> do/2. We will denote this set of vector 
fields by 5i. Let 5 = 5o U denote the family of tangential space vector fields and let T = 5 U {Dt} 
denote the family of space time tangential vector fields. 

Let the radial vector field be 

(6.6) R = ciy°'d/dy^, ci > 0 

Now, divR = n is not 0 but for our purposes it suffices that it is constant since what we need is that 
if divVF = 0 then ddvC rW = RdwW — W divi? = 0, where the Lie derivative Cr is defined in the 
next section. Let TZ = SU {R}. Note that IZ span the full tangent space of the space everywhere. Let 
U = {R} U {Dt} denote the family of all the vector fields construct above. Note also that the radial 

vector field commutes with the rotations; 

(6.7) [R,S]=0, SeSo 

Furthermore, the commutators of two vector fields in So is just ± another vector field in So- Therefore, 
for i = 0,1, let TZi = StU {i?}, % = StU {Dt} and Ui = SiU {i?} U {Dt}. 

Let U = {Ui}//.^ be some labeling of our family of vector fields. We will also use multindices 
I = (ii, ...,ir) of length |/| = r. so = Ui^ ■ ■ ■ Ui^ and C/j = Cu.^ • • • Sometimes we will write 

C/j, where U G So or I G So, meaning that G So for all of the indices in I. 

Note also that the vector fields U°‘{y)d/dy°' expressed in the x coordinates are given by U'^d/dx'' 
where [/* = U^'dx^/dy^. We here use the convention that indices refers to the components in 

the Lagrangian frame and indices i, ...,n refers to the components in the Eulerian frame. 
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7. Lie derivatives. 


Let us now introduce the Lie derivative of the vector field W with respect to the vector field T ; 

(7.1) = riL“ - {dcT‘^)W 

We will only deal with Lie derivatives with respect to the vector fields T constructed in the previous 
section. For those vector helds divT = 0 so 

(7.2) divlF = 0 ^ div£ThF = rdivlF-lFdivT = 0. 

The Lie derivative of a one form is defined by 

(7.3) CTaa = Taa +{daT‘^)ac. 

The Lie derivatives also commute with exterior differentiation, [Ct-, d] = 0 so if g is a function, 

(7.4) Crdaq = daTq. 

The Lie derivative of a two form is given by 

(7.5) CrPab = T(3ab + {daT^)Pcb + {dbT^)Pac. 

Furthermore if re is a one form and curlrcab = dwab = daW^ — d^Wa then since the Lie derivative 
commutes with exterior differentiation; 


(7.6) Tt CUrltCab = CUrlTytCab. 

We will also use that the Lie derivative satisfies Leibnitz rule, e.g. 

(7.7) CT{a,W^) = {CTa,)W^ + CAPacW^) = + PacCrW^. 


Furthermore, we will also treat Dt as if it were a Lie derivative and we will set 

(7.8) Cn, = A. 

Now of course this is not a space Lie derivative but rather could be interpreted as a space time Lie 
derivative in the domain [0,r] x Q. But the important thing is that it satisfies all the properties 
of the other Lie derivatives we are considering, such as divlF = 0 implies that divUtlF = 0 and 
Dtcuilw = cmlDtW, simply because it commutes with partial differentiation with respect to the y 
coordinates. The reason we use the notation (7.9) is that we will apply products of Lie derivatives and 

(7.9) and it is more efficient with the same notation. Furthermore 

(7.9) [Cd„Ct]=0 

this is because this quantity is and [Dt,T] = 0 for the vector fields we are considering, or it 

follows from (7.1) and that r“ = T“(y) is independent of t. 

17 



8. Commutators between Lie derivatives with respect to tangential 

VECTOR FIELDS AND THE NORMAL AND MULTIPLICATION OPERATORS. 

Note that the projection P defined in section 3 almost commutes with the Lie derivative with respect 
to tangential vector fields. In fact if denote the corresponding operator on one forms by P 

(8.1) PUa = Ua- daQ 

where q is as in (3.2) and Ua = QabU^■, then CrPua = — daTq. Since (7 = 0 on the boundary it 

follows that Tq = 0 there so the last term vanishes if we project again: 

(8.2) Ri^TPUa) = PPrUa 

We will need to calculate commutator between Lie derivatives with respect to tangential vector 
fields T and the operator Aj defined in section 3. Let Aj denote the corresponding operator taking a 
vector field to the one form 

(8.3) AfWa = gabAfW^ = -da{{dJ)W^ - q), 

Then since 

(8.4) Crda {{dcf)W^) = da {{daTf)W^) + da {{dJ)CTW-) 
it follows from ( 8 . 2 ) 

(8.5) PCrA^Wa = AfCrWa + A^^Wa 

Note that if / = p then it follows from (3.13) that the commutator is lower order. In fact p = 0 on 
the boundary implies that Tp = 0 on the boundary if T is a tangential vector field. Since VvP / 0 it 
follows Tp/p is a continuous function that is equal to VnTp/VnP on the boundary. Hence by (3.13) 

(8.6) \{W,AtpW)\ < \\VNTp/VNp\\L^i9n){W,AW) 

In view of (8.2) it follows that the multiplication operator Ma, defined by (3.16) in section 3, satisfies 
the commutator relation 

(8.7) PPtM^W = M^CtW + M^^^W, where M^Wa = PabM^W^ 
for a two form a. Let 

( 8 . 8 ) GT = MgT, g'^f, = CTgab, Ct = M^jT = Ctco 

We will also use special notation for the time derivatives of G: 

(8.9) G = Got = Mg, gab = DtPab 
and of A 

(8.10) At = Axp, A = At^ = A^tp 

In the following sections we will commute through products of vector fields = Ct^ ■ ■ ■ C,Ti where 
I = (L,..., v) and we will use the notation 

( 8 . 11 ) Ai = AtIp, Gi = Mgi, Gj = Moji 

where gL = Gr^tgab and = CipUJab- 
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9. Commutators between the linearized equation and Lie 

DERIVATIVES WITH RESPECT TO TANGENTIAL VECTOR FIELDS. 

We are now ready to commute tangential vector fields through the linearized equation and in the 
next section get the higher order energy estimates of tangential derivatives. Let T G T be a tangential 
vector fields and recall that [CT,Dt] = 0 and that if W are divergence free then so is CtW. Let us 
now apply Lie derivatives Cip = Cxi^ • • • where I = ...,ir) is a multi index, to the linearized 

equation (2.19) with an inhomogeneous divergence free term F vanishing to order r as t ^ 0: 

(9.1) gabW^ - da{{d,p)W^) = -daQ - {9ca " iO,a)W^ + 9abF\ = 0. 

which yields 

(9.2) {{d,TFp)cFw^) 

= -daT\ - 2 cIj^ ( 4 ^ {g,a - co,a))CFW‘^ + {CFg,^)CFF^ 
where we sum over all Ii +12 = I and 4/^ ~ introduce some new notation 

(9.3) Wi = CpW, Fi = CpF gif, = C^gab 4b = ^4ab, Pi = T^P, Qi = q 
and gif, = DtC^gab, Wj = DtWj etc. With this notation (9.2) becomes 

(M) chiJFv'L - chiM(8cn,)W'u) = -daq, - + Ufcs'ul 

Let us now project each term onto divergence free vector fields and also introduce some notation for 
the resulting operators 

(9.5) AiW^ = AxipW^, GiW^ = P{g^^gif,W^) 
and 

(9.6) GiW^ = P{g‘^^gif,W^), GjW^ = P{g^^ioif,W^) 

From now on we set 4^^ = 4/^ ''^hen I 2 ^ I and 4^^^ = 0 if /2 = /. Projecting each term onto 
divergence free vector fields we can now write (9.4) as 

(9.7) L^Wi = Wi + AWi + GWi - GWj = Pj - {Aj, Wi, + Gj, Wi, - GpWj, + Gj, Wj, + Gj,Pi,) 
Here Gj, Gj and Gj are all bounded operators. By (3.16)-(3.17): 

(9.8) \\GjW\\ < \\CU\L^(n)\\W\\, \\GjW\\ < ||£4|U»(f,)||W|| 

The terms Gi, Wj, are easy to take care of by also including time derivatives up to highest order in our 
estimates since I/ 2 I < |7| — 1. AWi itself will be included in the higher order energy, which is just going 
to be a sum of terms of the form (4.4) with W replaced by Wi for |/| < r. However, we also have to 
deal with ApWj, since Aj, is an operator of order 1. Since I/ 2 I < |/| — 1 < r — 1 in the terms ApWj, 
and since the energy will give us Wi for all |/| < r we in particular will have an estimate for Wi, which, 
using the equation (9.7), up to terms of lower order is —AWj,. Since Aj = A^jp it follows from (3.13) 
that 

(9.9) \{U.AjW)\ < \\VNT-^p/VNP\\L^idii){U,AU)^/^{W,AW)P\ 

However this does not imply that the norm of Aj is bounded by the norm of A. Therefore we have to 
deal with these terms with Aj, in an indirect way, by including them in the energy and using (9.9). 
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10. The a priori energy bounds for tangential derivatives. 

To obtain estimates for higher derivatives we apply tangential vector fields to the equation and get 
similar equations for higher derivatives. However, there are a some commutators coming up that we 
have to deal with. One can be dealt with by adding a lower order term to the energy and another 
commutator one deals with by also considering higher time derivatives. The main point is however that 
commutators with the normal operator can be controlled by the normal operator through (9.9). Let 
Wt = CtW, Ft = CtF and let Gt-, Ct and At be as in (8.8) and (8.10). By (9.7) 

( 10 . 1 ) LiWt = Ft- AtW - GtW + GtW + GtW + GtF 

The terms one has to deal with are AtW and GtW. Let E = E{W), where E{W) is given by (4.4), 

(10.2) Et = E{Wt) = {Wt,Wt) + {Wt,{A + I)Wt), and Dt = 2{Wt, AtW) . 

Dt is lower order compared to Et since by (9.9) it is bounded by a constant times y/ET^/E and we 
already have an estimate for E in (4.10). We will add Dt to the energy Et to pick up the commutator 
At between Ct and A. By (4.8) 

(10.3) Et + Dt = 2{Wt, LiWt) + 2{Wt, Wt) - {Wt, GWt) + {Wt, AWt) + {Wt, GWt) 

+ 2{Wt, AtW) + 2{Wt, AtW) + 2{Wt,AtW) 

= 2(1Tt, AtW) + 2{Wt, GtW) + 2{Wt,E + GtE) 

+ 2{Wt, -GtW + GtW + Wt) - {Wt, GWt) + {Wt, AWt) + {Wt, GWt) + 2{Wt,AtW) 


Here, the terms on the last row are bounded by Et and E using (9.8) and (9.9). The only terms that 
remains to control are 2(Wt,GtW) and 2{Wt,AtW). These terms are controlled by simultaneously 
consider one more time derivative, i.e. if T = Dt, and estimate energies for these. 

Let us now define higher order energies. Let 
(10.4) El = E{Wi) = {Wi, Wi) + {Wi, {A + I)Wi), Wi = C^rW. 

With notation as in the previous section we have by (4.8) and (9.7) 


(10.5) El = 2{Wi,DtWi + AWi + GWi - GWi) 

+ 2{Wi,Wi) - {Wi,GWi) + {Wi,AWi) + {Wi,GWi) 


= -2c 


hh 


Wj, Aj^ Wi,) + {Wi, Gi, Wi,)-{Wi, Gi, Wi,) + {Wi, Gi, Wi,) + {Wi, Gi,Ej, 

+ 2{Wi,Ej) + 2(1T/, Wi) - {Wi,GWi) + {Wi,AWi) + {Wi, GWi) 


To deal with the term {Wj, Aj^Wj^) we introduce 

(10.6) Di = 2cY^-{Wi,Ai^Wi,) 

Then 


(^^i,AtWi,) + {Wi,AtWi,) + {Wi,AtWi,)) 
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(10.7) 


Dj = 



and hence 


(10.8) Ej + Dj — 

+ 2{Wi,Fj) + 2{Wi, Wi) - {Wi,GWi) + {Wi,AWi) + {Wi, GWi) 

We have hence replaced the bad term by two terms that we can control by (9.9). Furthermore, we can 
also bound Dj itself using (9.9). 

For a two form a and a function q vanishing on the boundary let 

(10.9) ||a||oo = II |a| ||L<»(f7), l|99lloo,p-i = W'^Nq/'^NPho-idn) < II^^IU/cq, 
and for a vector fields W let 

(10.10) {W)a = {W,AW)^/‘^, \\W\\ = 

With this notation it now follows from (10.8) and (9.8)-(9.9) that 

(10.11) Ej + Dj < 2{Wj)Acj^^^ {\\dph lloo,p-i {Wi^)a + \\dpii\\oo,p-^{Wi^) a^ 

+ 2|| Wp|| (dl/^ lU + ||c^'^ \\o.)\\Wi, II + 11/^ lUII^/, II + 11/^ \U\Fj, ll) 

+ ||lFp||(2||Fp|| + 2||Wp|| + ll^lloollW/ll) + llf^/ll IliilUllf^/ll + {Wj)A\\dp\U,,-.{Wi)A 

Furthermore 

(10.12) \Dj\<2{Wi)Ac}^^^\\dpi,\U^p-i{Wj,)A 



If I £ F and |/| = r then with the notation in Definition 10.1 we obtain from (10.11) and (10.12); 


(10.16) 


\El + i)z| < GE-^ Y. i^r-s + ll^ll^-s), 


\Dj\<GE^Y^^^r-l- 
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If we integrate the first inequality from 0 to t using that Ei(0) = Di{0) = 0 and the second inequality 
we get with a constant depending on = supo<,-<T nj (r) 

(10.17) Ei<CEJeJ_^ + C ! eJ{eJ + \\F\\'^)dT 

Jo 

If we sum over \I\ < r and divide by Er{t) = supQ<.^<j EJ (t) we get for some other constant 

(10.18) Er<CEr-i + C f [Er + \\F\\J) dr. 

Jo 

Hence with Mr{t) = E^ dr, we get 

(10.19) -CMr <CEr-i + C [ \\F\\J dr 

at Jq 

Multiplying by the integrating factor and integrating from 0 to t we see that is bounded by 
some constant depending on t < T times the right hand side and hence it follows that for some other 
constant 


( 10 . 20 ) 


Er “E CEy—i 



F\\J dr 


Since we already proved a bound for Eq in (4.10) it inductively follows that: 

Lemma 10.1. Suppose that x,p G C'’'+^([0, T] x H), p\q^ = 0, < —cq < 0 and divH = 0, 

where V = Dtx. Suppose that W is a solution of (9.1) where F is divergence free and vanishing to 
order r as t ^ D. Let EJ be defined by (10.14). Then there is a constant C depending only on the 
norm of {x,p), a lower bound for cq and an upper bound for T, such that if Ej (0) = 0, for s <r, then 

(10.21) eJ {t)<C [ \\F\\J dr, for 0<t<T 

Jo 


11. Estimates of derivatives of a vector field in terms 

OF THE CURL, THE DIVERGENCE AND TANGENTIAL DERIVATIVES. 

In this section we show that derivatives of vector fields can be estimated by derivatives of the curl, 
the divergence and tangential derivatives. First we prove the basic estimate in the Euclidean coordinates 
in Lemma 11.1 below. This estimate it is not invariant and so in Lemma 11.2 we express it in terms of 
Lie derivatives which is invariant. 

Lemma 11.1. IFe have 


(11.1) |5a| < Cnd curla| + I diva| + Xlsgsl'S'al), cmlaij = diOj — djOi 


diva = d'^^diUj 


for a one form ai in the Eulerian frame, where Cn only depends on the dimension n. Here the norms 
are the Euclidean norms, \da\ = \diaj p. 
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Proof of Lemma 11.1.. Since S span the full tangent space in the interior when the distance to the 
boundary d{y) > do we may assume that d{y) < do. Let 17“ = {y, d{y) > a} and let Vf^ be the image 
of this set under mapping y x{t, y). Let N the exterior unit normal to dDf'. Then — N'^N^ 

is the inverse of the tangential metric. Since the tangential vector fields span the tangent space of the 
level sets of the distance function we have S’^S^OiOj, where here S'* = S'“9x*/(9y“. 

We claim that for any two tensor fdij-. 

(11.2) + |/3p + (tr/3)2) 

where fdij = jdij — fdji is the antisymmetric part and tr/3 = (3ij is the trace. To prove (11.2) we 
may assume that (3 is symmetric and traceless. Writing 5*1 = g*l + A^* W we see that the estimate for 
such tensors follows from the estimate N^N^PkiPij = < nq^^PkiPij■ 

(This inequality just says that {tic{QP)p < nti{QPQP) which is obvious if one writes it out and use 
the symmetry. ) □ 

The inequality (11.1) is not invariant under changes of coordinates so we want to replace it by an 
inequality that is, so we can get an inequality that holds also in the Lagrangian frame. After that we 
want to derive higher order versions of it as well. The divergence and the curl are invariant but the other 
terms are not. There are two ways to make these terms invariant. One is to replace the differentiation 
by covariant differentiation and the other is to replace it by Lie derivatives with respect to the our 
family of vector fields in section 6. Both ways will result in a lower order term just involving the norm 
of the one form itself multiplied by a constant which depends on two derivatives of the coordinates. 

Definition 11.1. Let ci be a constant such that 

(11.3) {\gab\ + < cj, Idx/dyl"^ + \dy/dx\‘^ < cj 

a,b 

and let Ki denote a continuous function of ci. 

We note that the bound for the Jacobian of the coordinate and its inverse follows from the bound 
for the metric and its inverse and the bound for the Jacobian and its inverse implies an equivalent 
bound for the metric and its inverse with cf multiplied by n. All our constants in what follows in this 
section will depend on a bound for ci and we will denote such a constants by Ki. 

Lemma 11.2. In the Lagrangian frame we have, with = gabW^, 


(11.4) 

\CuW\ < iLi(|curlW| 

+ \d\yW\ + Y:ses\^sW\ + [g]pW\), 

u en, 

(11.5) 

\CuW\ < ATi^lcurl W| 

+ |divW| + + [g]i\W\), 

U eU 

where [^Ji 

= 1 + \dg\. Furthermore 




(11.6) \dW\ < K, (\CrW\ + + |11^|) 


When d{y) < do we may replace the sums over S by the sums over So and the sum over T by the sum 
over Tq. 

Proof of Lemma 11.2. (11.5) follows directly from (11.4) by adding the time derivative to the right hand 
side. We will show that (11.4) in the Eulerian frame follows from (11.1) and then it follows directly 
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that (11.4) holds in Lagrangian frame as well since everything is invariant. Let Z'^ = Then 

LijZ'^ = [/Z* — {dkU^)Z^, where W = U°'dx'^/dy°‘ are the components of the vector field U expressed 
in the Eulerian frame. Now transforming to the Lagrangian frame, partial differentiation becomes 
covariant differentiation. {dkU^){dx^/dy°'){dy^/dx^) = VaU^, where VaU^ = daU^ + and 

= 9‘'‘^ (daObd + dbOad “ ddgab) /2 = {dy^/dx^)dadbx'' are the Christoffel symbols. Since \daU^\ < C 
it follows that < C\g\i. That we may replace S by Sq close to the boundary follows from the proof 

of Lemma 11.1. (11.6) follows since TZ span the tangent space and |£[/1T“ — UW°‘\ = \{dcU°')W‘^\ < 
C\W\. □ 

We are now going to derive higher order versions of the inequality in Lemma 11.2. We want to 
apply the lemma to W replaced by CjjW. Then in our applications the divergence term vanishes and 
as we shall see later on we will be able to control the curl of {C^fW^a = ^uidabW^) which however is 
not the same as the curl of {^LyW)a = gabZZfjW^ but the difference is lower order and can be easily 
estimated. Let us hrst introduce some notation: 

Definition 11.2. Let /? be a function, a one or two form or vector field, let V be any of our families of 
vector fields and set 

(11,7) Y. i/Jib-'i/Jii.. i^i5’ = i, 

|J|<S, JeV Sl+...+Sfc<M,Si>l 


In particular |/3|^ is equivalent to \dyfi\ and |/3|(^ is equivalent to Yl\a\+k<r \^tdyfi\- 

Lemma 11.3. With the convention that |curlIT|Yi = |divIT|Yi = 0 rce have 

r 

(11.8) < iLi (I curlWi:".! + | divW|:"_i + | \g\f\W\^..s), 

S = 1 

r 

(11.9) |W|^ < E (|curllT|:^_i_, + |divW|:^_i_, + |W|f_,). 

s=0 


The same inequalities also holds with TZ replaced by U everywhere and S replaced by T: 

r 

(11.10) \W\^ < iLi( |curlW|"_i + |divW|"_i + \W\J + \g\^\W\^_,,), 

S = 1 

r 

(11.11) \W\^ < E ^9]s {I curlW|(^_i_, + I divW|(^_i_, + \W\J_,). 

s=0 


Proof of Lemma 11.3. We will first prove (11.8) We claim that 

(11,12) E \C{jW\<Ki (IcurlT^I + IdivT^ITl + [5 ]i|£^IT|)+ iLi E Ts"'\ 

\I\=r,Uen \J\=r-l,UeTl \l\=r,ses 

First we note that there is noting to prove if d{y) > do since then S span the full tangent space. 
Therefore, it suffices to prove (11.12) when d{y) < do and with S replaced by So and TZ replaced by 

24 



TZq. Then (11.12) follows from (11.4) if r = 1 and assuming that its true for r replaced by r—1 we will 
prove that it holds for r. If we apply (11.4) to where | J| = r —1, we get 

(11.13) \CuCijW\ < i^i(|curl£^lT| + |div£^ir| + ^ \CsCiW\ + [g]^\CijW\). 

ses 

If C(j consist of all tangential derivatives then it follows that |£f/T^lT| is bounded by the right hand 
side of (11.12). If does not consist of only tangential derivatives then, since [£^,£ 5 ] = £[_r, 5 ] = 0, 
if S' G 5o, we can write CsC'^W = Cs'W, for some S' G Sq. If we now apply ( 11 . 12 ) with r replaced 
by r — 1 to £ 5 'IT, (11.12) follows also for r. 

In (11.8) we have C\, curllT = curl£y^lT which however is different from curl C\jW. We have; 

(11.14) £^^1T„ = CUga.W'') = -ga,C-ijW^ + where g^ = C^g^, 

where the sum is over all J\ + J 2 = J and Cj^j^ = 1 for | J 2 I < | J| = 0 if J 2 = J. It follows that 

( 11 . 15 ) |curl£^IT-curl£^IT| + VU < |J|, 

where the partial derivative can be estimated by Lie derivatives. (11.9) follows by induction from (11.8). 
Finally, (11.10) follows from (11.12) and (11.15). In fact, applying (11.12) to IT replaced by C^d^W we 
see that (11.12) holds also for TZ replaced by U and S replaced by T and (11.15) also holds for U ^U. 

12. The estimates for the curl and the normal derivatives. 

Note that in section 10 we only had bounds for the derivatives that are tangential at the boundary, 
as well as all derivatives in the interior since S span the full tangent space in the interior. We will now 
use estimates for the curl together with the estimates for the tangential derivatives to get estimates also 
for normal derivatives close to the boundary. Let 

(12.1) Wa = gabW^, and curl Wab = daWb - dbWa- 

Then we have 

( 12 . 2 ) Dt{gabW>') - da{{d,p)W^) = -daq + COabW^ + F^ 

Note that (12.2) can also be formulated as 

(12.3) Dtw + AW - CW = F 

where the underline as before means that we lowered the indices so the result is a one form. Note here 
that w is not equal DtW so the notation is slightly confusing. But what we mean is that we think of 
IT as a vector field and take the time derivative as a vector field which results in IT and then w is the 
corresponding one form obtained by lowering the indices. We obtain 

(12.4) A curl Ab = (5e LVab)W^ - co.bdaW'' + uj^adbW^ + 

Since DfWa = gabW^ + gabW^ and dagbc - dbgac = dcOJab we also obtain 

(12.5) A curl Wab = curl Wab + idc(^ab)W‘' + gbcdaW" - gacdbW. 
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Since di\W = di\W = 0 it follows from Lemma 11.2 and (12.4)-(12.5) that 

(12.6) |L)tCurlzc| < iLi |t(;| (| curltc| + + b]il^l) + l^wHVLl + | curlFj 

(12.7) lAcurltcl < Icurltcl + iLil^l (| curlu;| + Mil^l) + |5a;||lL| 

Since we already have control of the tangential derivatives S by section 11 this obviously gives us 
control of curltc and curlrc as well and once we have control of these we in fact control all components 
by Lemma 11.3 again. The norms will be measured in since we have control of the norms of the 
tangential components. We will now derive higher order versions of the inequalities (12.6)-(12.7) using 
the higher order version of Lemma 11.2, i.e. (11.9) in Lemma 11.3. 

We must now get equations for the curl of higher derivatives as well. Applying to (12.4)-(12.5) 
gives, since the Lie derivative commutes with the curl, 

( 12 . 8 ) DtCuiWljWab = + + {cvj\CijF)ab, 

where = C^j lv and 

(12.9) A cnilC^Wab = curlT^Ab + + cj,j, 

where = C(jDtgab- Let us make a definition: 

Definition 12.1. Let /3 be a two form. With notation as in Definition 11.2 we set 

(12.10) {Wi\)l = E Il'li’i'Sll- 

s+r<w 


Using Lemma 11.3 and Lemma 11.2 it follows that: 

Lemma 12.1. With notation as in Definition 11.1 and Definition 12.1 and the convention that 
I curllUl^j = |divlU|):(;^ = 0 we have 

r 

(12.11) \DtCuvlw\'^_i < Ki^2 (b]|<^l)J^s(|curlu;|^;^ + |divlU|^i + |lU|f) + |curlF|^;^ 

s=0 

r 

(12.12) |Acurlu;|^_i < A X] (I curlu;|f_i + |divlU|^_i + |W|f) + |curlu;|^_i 

s=0 

The same inequalities hold with TZ replaced by U and S replaced by T. 

Proof of Lemma 12.1. Let us hrst prove (12.11). The first terms in the right hand side of (12.8) are by 
Lemma 11.3 bounded by a constant times 

r r u 

(12.13) < A'. + |divivif_,+|ivif) 

u=0 u=0s=0 

The proof of (12.12) uses the same argument and that that dcWab = daPbc — dhPac- D 
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Let us now introduce some new norms and some new notation: 

Definition 12.2. For V any of our families of vector fields let 

(12.14) \\wf^ = \\wmv^(n)= (/ \ChW{t,y)fKdyy^\ 

\l\<r,lev 

and 

(12.15) Cy = ^ f I curl£^rcp + I curl£^rcp , Cq =0. 

\J\=<r-i,Jev 

Note that ||lF(t)|| 7 ^>-(o) is equivalent to the usual Sobolev norm in the Lagrangian coordinates. 

Definition 12.3. For V any of our families of vector fields and for fi a function, a 1-form, a 2-form or a 
vector field let \P\^ be as in Definition 11.1 and set 

( 12 . 16 ) Ili3||^^ = |||/3irilt-(m. |[9lC= E NiLco-■ • [[9K== = 1 

Si + ...+Sk<S, Si>l 

where the sum is over all combinations with > 1. Furthermore, let 

(12.17) rn^ = M%oo^ '^r = X] (11^1+ Il“.“)' 

sH-w<r 


Let = llcurlF It now follows from Lemma 12.1 that 


(12.18) 


K^Y,rh^-s{C^ + El)+Fi: 


where EJ is the energy of the tangential derivatives defined in section 10. Hence 


t , r-l 


(12.19) 


< i^ielo K^rh.dr / ^ ^ + F^^ ) dr 


Since we already proved a bound for EJ in Lemma 10.1 it inductively follows that C\! is bounded. 
Note that, if r = 1 the interpretation of (12.19) is that the first sum is not there. By Lemma 11.3: 


( 12 . 20 ) 


W{t)\\u HO) + \\W{t)\\u -(O) < Ki 


Hence we have: 


Lemma 12.2. Suppose that x,p G C'’'+^([0, T] x H), = 0, < —cq < 0 and divF = 0, 

where V = Dtx. Then there is a constant C = C{x,p) depending only on the norm of {x,p), a lower 
bound for cq and an upper bound for T, such that if EJ (0) = (0) = 0, for s < r, then 


( 12 . 21 ) 


U , tpT 


+ Ef<C \\E\\‘f dr, for 0<t<T. 
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13. The smoothed out normal operator. 


In order to prove existence we first have to replace the normal operator ^ by a sequence A'^ of 
bounded symmetric and positive operators that convergence to A, as e —> 0. The boundedness is 
needed for the existence and the symmetry and positivity is needed to get a positive term in the energy. 
Furthermore the commutators with Lie derivatives with respect to tangential vector fields as well as the 
curl have to be well behaved. Let p = p{d) be a smooth function of d = d{y) = dist (y, 911), such that 

(13.1) p' > 0, p{d) = d for d < 1/4 and p{d) = 1/2 for d > 3/4. 

Let x{p) be a smooth function such that 

(13.2) X^(/o) > 0, x(d) = 0) when p < 1/4, and y;(yo) = 1, when p>S/4: 

For a function / vanishing on the boundary we define 

(13.3) = P{- g‘^^Xe{p)db{fp-\d,p)W^)) 
where Xeip) = x{p/^)- Then if we integrate by parts we get 

(13.4) {U,A}W)= [ fp-\',{p)iU‘'dap){W’’dbp)dy 

Jn 

from which it follows that A’y is symmetric and 

(13.5) A^ = A^p>0, i.e. (1F,T^1T)>0, if p>0 
It also follows that another expression for is 

(13.6) A}1T“ = P{g‘^^x'e{p)idbP)fp-\d,p)W^) 

A'^ is now for each e > 0 a bounded operator 

(13.7) \\A^W\\<C\\VNp\\L^e-^\\W\\ 

since Xe — VP~^\dp\ < C\\VNp\\L°°{dn)- Iii general, since the projection is continuous on 

id’'(ll), see (3.6) and (3.8), if the metric and pressure are sufficiently regular we get 

k k 

(13.8) ^\\DiA^W\\Hr(^n)<C,,r,k^\\DlW\U^n). 

j=0 j=0 

Moreover 

(13.9) A^U ^ AU, in if U £ H^{n) 

In fact, the projection is continuous in the norm and XeF ^ T in if T G L^. It follows that 

(13.10) P{g‘^\,{p)db{pp-\d,p)U‘^)) ^ P{g<^^db{pp-\d,p)U‘^)) = P{g<^^db{{d,p)U‘^)) 


since p p ^dcP = dcP on the boundary. 
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We will now calculate the commutators with the Lie derivative Ct with respect to tangential vector 
fields T. As before the inequality 

(13.11) \{U,A%W)\ < 

hold, where = {y G 11; d{y) > e}. In fact, it suffices to take the supremum over the set where 
d{y) < £ since x'e = 0) when d{y) > e. The only difference with (3.13) is that now the supremum over a 
small neighborhood of the boundary instead of on the boundary. The positivity properties (13.5) and 

(13.11) for will play the role that (3.12) and (3.13) did for A. In particular, since p vanishes on the 
boundary, p > 0 in the interior and VtvP < —cq < 0 on the boundary it follows that p = Dtp vanishes 
on the boundary and p/p is a smooth function. Therefore 

(13.12) A^ = A^. satishes \{W,A^W)\<\\p/p\\L^^^^\^e){W,A^W) 

Here A^ is the time derivative of the operator A^, considered as an operator with values in the one 
forms. It will show up in the energy estimate for the e smoothed out equation in the next section. 

The commutators between Aj and Lie derivatives with respect to tangential vector fields are basi¬ 
cally the same as for A. Note that 

(13.13) r(i = o, if rGro = 5oU{A} 
where Sq are the rotations. Hence if T is any of these vector fields we have 

(13.14) P{g^<^CT{gabA}W^)) = A)CtW^ + A^^fW^, 

However, in order to get additional regularity in the interior we include the vector fields 5i that span 
the tangent space in the interior. The vector fields in 5i satisfy 

(13.15) Sp = CsP = 0, when d<do/2 

Since x'eip) — 0 when d > £ the commutator relation (13.14) above is true for these as well if we assume 
that £ < do/2. 

It remains to estimate the curl of A^. Whereas, the curl of A vanishes this is not the case for the 
curl of A^. It will however vanish away from the boundary. With A^Wa = PabA^W^ we have 

(13.16) A^Wa = -Xe{p)da{p P~^{dcP)W'') - daPl 

for some function qi vanishing on the boundary and determined so the divergence vanishes. Since the 
curl of the gradient vanishes and x'e{p) = 0 when d > £ we have 

(13.17) curl A'^ITab = 0, when d{y) > £ 

14. The smoothed out equation and existence of weak solutions. 

The £ smoothed out linear equation; 

IT“ + A^IT“ + GITe-CHT = F“, =0 

t=o t=o 
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(14.1) 



is just an ordinary differential equation for [We,We) on the space of divergence free vector fields in 
since all operator are bounded so existence follows in L^(f7). In fact its an ordinary differential 
equation in the Sobolev spaces by (13.8). To get additional regularity in time as well we apply 

more time derivatives using (13.8) and (3.8) and that the initial conditions for these vanishes as well 
since we constructed F in (14.1) so it vanishes to any given order. If initial data, encoded in F, are 
smooth, we hence have a smooth solution of the e approximate linear equation. 

Now we want to use the existence and estimates for the e smoothed out linear equation and pass 
to the limit as e —> 0 to get existence for the linearized equation. Will show that W weakly in 

L^, where W G for some large r. From the weak convergence it will follow that IF is a weak 

solution and then from the additional regularity of IF it will follow that in fact its a classical solution 
and hence that the a priori bounds in the earlier section hold. 

Of course the norm of tends to infinity as e ^ 0 but since it is a positive operator it can be 
included in the energy. The energy will be the same as before with A replaced by A^, so (4.4) becomes 

(14.2) = {We,We) + {We,iA^ +I)We) 


The time derivative of the first term is the same as (4.5) with IF replaced by IF^. Since Dt d = 0 it 
follows from taking the time derivative of (13.4), with f = p, that 

(14.3) j^{We,A^We) = 2{We,A^We) + (IF,, A^IF,), 
where the last term is bounded by (13.12). Hence by (4.7)-(4.9): 

(14.4) \E^\ < (l + + \\{Dtp)/p\\L^^n))E^ + 2V^\\F\\ 


from which we get a uniform bound for 0 <t <T independent of e: E’^{t) < C. 

Since ||II4|| < C we can now choose a subsequence IF,^ ^ IF weakly in the inner product. We 
will show below that the limit IF is a weak solution if the equation. Multiplying the e smoothed out 
equation by a smooth divergence free vector field U that vanishes for t > T and integrating by parts 
we get 

rp rp 

(14.5) [ [ 5ab(i/“ + H^t/“ + G17“-C'f7“-Ct/“)lFj’d2/(it= [ [ gahU^F'^dydt 

Jo Jn ^ ^ Jo Jn 

where C'lF'^ = P(^g°'‘^(jjcbW^), since A^ and D^ + BDt are symmetric and the adjoint of CDt is CDt — C. 
We proved in the previous section that A^U converges to AU strongly in the norm if U is in . Since 
IF,^ ^ IF weakly this proves that we have a weak solution IF of the equation: 

T T 

(14.6) / [ gabiu'" + AU^^ + GiF^ - Cil^^ - CU‘^)w^ dydt = [ [ gabU'^ F^ dydt 

Jo Jn ^ ^ Jo Jn 

for any divergence free smooth vector field U that vanishes for t >T. Furthermore since IF, is divergence 
free, we have 

T 

(14.7) / [ idaq)W^ dydt = 0 

Jo Jn 

for any smooth q that vanishes on the boundary and hence 


(14.8) 


{daq)W<^ dydt = 0 


so IF is weakly divergence free. 
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15. Existence of smooth solutions for the linearized equation. 

Now that we have existence of a weak solution we will prove that we have additional regularity and 
in fact that, W,W€ //’’(fl) for any r > 0. It then follows that we can integrate by parts again in the 
above integrals and conclude that 

T 

(15.1) [ [ qdaW‘^dydt = 0 

Jo Jn 

for any smooth function q that vanishes on the boundary. Hence W is divergence free. Furthermore 

rp rj-i 

(15.2) / / gabU‘^(w^ + AW^ + GW^-CW^)dydt= [ [ gabU^F^dydt 

Jo Jn ^ ^ Jo Jn 

for any smooth divergence free vector field U that vanishes for t >T. But in fact since W is divergence 
free it follows that + AW^ + GW^ — CW^ is divergence and since by construction F is divergence 
free as well it follows that (15.2) holds for any smooth vector field U that vanishes for t >T. We then 
conclude that 

(15.3) + AW^ + GW^ - CW^ = F\ divIF = 0 

It therefore only remains to show that W G H’’(H). We must show that we have uniform bounds 
for the £ smooth out equation similar to the a priori bounds for the linearized equation. 

The uniform tangential bounds for the e smoothed out equation follows the proof of the a priori 
tangential bounds in section 10. The proof is just a change of notation. Let 

(15.4) El = {Wei, Wei) + {Wei, (H + I)Wel), Wei = C^rWe- 

If £ < do then the commutator relation for A^, (13.14), is exactly the same as for A, (8.5). Furthermore 
the positivity property for A^ only differs from the one for Af hy that the supremum over the boundary 
in (3.13) is replaced by the supremum over a neighborhood of the boundary where d{y) < e in (13.11). 
Hence all the calculations and inequalities in sections 10 and 12 hold with A replaced by A’^, if we replace 
the supremum of V]v<?/VjvP over the boundary in (10.9) by the supremum of q/p over the domain 
where is given by (15.6). Therefore we will arrive at the energy bound (10.21) for EJ replaced by 

(15.5) ^ 

\i\<r,ier 

i.e. Lemma 10.1 hold for EJ replaced by EJ'^ with a constant independent of e. Note that, this is 
where we need to have vanishing initial conditions and an inhomogeneous term that vanishes to high 
order when t = 0 so that also the higher order time derivatives of the solution of (14.1) vanished when 
t = 0. If the initial conditions for higher order time derivatives were to be obtained from the e smoothed 
out equation, then they would depend on e and so we would not have been able to get a uniform bound 
for the energy, EJ' ^. 

The bound for curl is very simple since by (13.17) the curl of Ae vanishes in 


(15.6) 


= {y; dist(y,9H) > e}, 
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it follows that all the formulas in section 12 hold when d{y) > £. This follows from replacing A in (12.3) 
by and using that the curl of this vanishes for d{y) > e. Since all the estimates used from section 11 
are point wise estimates we conclude that (12.11)-(12.12) hold for W replaced by when d{y) > e. 
Let 


(15.7) 


= 


E 

I Jeu 




curlT^tCep + I curlT^tCgp 


1/2 


With replaced by and EJ replaced by we get exactly the same inequalities as before 

(12.18)-(12.19), since these were derived from the point wise bounds in section 11. Furthermore, the 
inequality (12.20) hold as well if we replace the norms by 


(15.8) 


\\wm 




\l\<r,leu 




1/2 


Therefore we conclude that the inequality in Lemma 12.2 hold with a constant C independent of e if 
we replace the norms by (15.8): 

Lemma 15.1. Suppose that x,p G C'’'+^([0, T] x fl), p\q^ = 0, '^np\q^ < —Co < 0 <^nd divF = 0, 
where V = Dtx. Suppose that is a solution of (If-l) where F is divergenee free and vanishing 
to order r as t ^ 0. Let EJ'^ he defined by (15.5). Then there is a eonstant C depending only on 
the norm of {x,p), a lower hound for cq and an upper bound for T, but independent of e, sueh that if 
-^^’^(0) = = 0, for s <r, then 

(15.9) \\W,{t)\\u^^ne) + mmu^{n^)+E^'%t)<C f\\F\\)^dT, for 0<t<T. 

Jo 


It therefore follows that the limit W satishes the same bound with replaced by Li, and so the 
weak solution in section 14 is in fact a smooth solution. 

16. The energy estimate revisited and the proof of the theorem. 

In section 10 we estimated the energies of the tangential derivatives without using the estimate 
of the normal derivatives coming from the curl. This was necessary to get uniform bounds for the e 
smoothed out equation since in that case we could not estimate the curl close to the boundary. The 
drawback was that instead we had to include all time derivatives as well in the energy. However, now 
that we have existence we can obtain other bounds for the linearized equation directly. In section 
9 we calculated the commutator between the linearized operator, considered as an operator from the 
divergence free vector fields to the one forms, and Lie derivatives with respect to tangential vector fields, 
and then projected the result back onto the divergence free vector fields. This was needed because the 
commutator between Lie derivatives and the operator A considered as an operator with values in the 
one forms is better behaved. However, the drawback is that the commutator with the second time 
derivative, considered as an operator with values in the one forms, involves second time derivatives, 
which is why we had to include all the time derivatives. Now we will instead commute through directly 
with the operator from the divergence free vector fields to the divergence free vector helds. Let us then 
also consider the original setting with non vanishing initial conditions and an inhomogeneous term; 

IT“ - g'^’’db{{d,p)W^ - q,) = -g^^iigcb - 
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where qi and q2 vanishes on the boundary and are chosen so that each term is divergence free. The 
second term on the left is AW^ and the term in the right is — + CW^'. Let us now first calculate 

the commutators with A and tangential vector fields. 

(16.2) Cs[g‘^^db{{d,p)W--qi)') 

= iCsg‘^'’)db{id,p)W^ - gi) + g‘^^dk{id,Sp)W^ + {d,p){CsWr - Sqr) 

where C-sg^"^ = —g°'‘^g’^^gfd, gf^ = Csgcd- Projecting each term onto divergence free vector fields: 

(16.3) CsAW =-GsAW + AsW + ACsW, 

where As = Asp and Gs = MgS is the operator GslT“ = P{g^^^g^b^’’) ■ Expressed differently 

(16.4) [Cs,A]W = {As-GsA)W 

Although Gs is a bounded operator, all the positivity properties of A are lost and the best we can say 
is that GsA is an operator of order 1. The operator As is also of order 1 but in section 10 we used 
the positivity property to estimate it in terms of A which we controlled by the energy. It remains to 
calculate the commutator with Gs and C, which basically are the same. 

(16.5) CtGsW^ = Cr^g'^^glW- - dbq)) 

= {Crg'^^) {glw'^ - dbq) + g'^\CTgl)W‘^ + g'^^glCrW^ - g'^^dbTq 

Projecting each term onto the divergence free vector fields we arrive at 

(16.6) [Ct,Gs]W = {Gts — GtGs)W, 

where GtsW^ = P{g‘"^glfW‘') and gff = CxCsgbc- 

In general using (16.4) and (16.6) to commute through we get for some constants 

(16.7) C^AW - AC^sW = • • • Gi,_,Aj,_^Wi, 

where the sum is over all combinations with Ii + ... + Ik = I, with k > 2, and \Ik\ < |7|. Here 
= MgjlT“ = P{g‘^^gcbW^), where g;^^ = Cggaci Aj = Agjp and ITj = Similarly we get 

the commutators with G and G 

(16.8) C^sGW - GC^sW = • • • Gj,_,Gi,_,Wj, 

(16.9) C^sGW - GC^W = ■ • • Gi,_,Gi,_,Wi, 

The only thing that matters is that these are bounded operators, and in fact they are lower order since 
\Ik\ < |/|. Hence we obtain 

LiW = Wi + AWi + GWi - GWi = Hi 
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where 


(16.11) Hj = Fj+ 

+ • • Gi,_,Gi,_,Wi, + • • Gi,_,Gi,_,Wi, 

where \Ik\ < |/| in the right hand side. Here Fi = C^F. As before, let 

(16.12) Ei = {Wi,Wi) + {Wi,{A + I)Wi) 

where we now only consider Wi = C^gW with S £ S. The energy estimate is like before and we only 
have to be able to estimate the norm of the right hand side of (16.10). The terms on the second 
row of (16.11) are obviously bounded by Ej for some |J| < |/|. In fact they are even lower order 
since we have strict inequality. Therefore it only remains to estimate the term on the right in the first 
row. |Ifc| < |/| but A/j. is order one and it contains derivatives in any direction so that term has to be 
estimated by the and so it does not directly help to have an estimate for 

for all tangential derivatives S. However the estimate of the tangential derivatives together with the 
estimates for curl in Lemma 12.1 gives the required estimate. 

Let G^ be defined (12.15), let be defined by (10.13) and let and m’^ be as in Definition 
12.3. Then by Lemma 11.3 we get the inequality corresponding to (12.20); 

r 

(16.13) ||iy||^ + llli^ll^ < K,J2^f-s{C^ + Ef), where ||1T||, = \\W{t)\\n^^n) 

s=0 

Since the projection has norm 1, ||GjlT|| < || 5 ''^||oo||lT||. It follows that 

(16.14) ||G,, • • • Gj,_,Gj,_,WiJ < ll/Mloo • • • 11/'=-^Hoc < rhf_J\W\U 

(16.15) ||G,, • ■■Gi,_,Gi,_,WiJ < ll/Mloo • • • 11/'=-^Hoc ||/''-Mloo \\WiJ < m:"_,||lT|U 

where s = \Ik\ < r and r = |/|. Let 

(16.16) Pf = [b]]?-^,oo X] \\dS-^p\\L--{dn) 

s=0 |J|<s+l,Je5 


Since Aj = Agjp it follows form (3.15) that 

(16.17) ||G,,. • < ll/MU- • • ||5'''-^||oo 

By (4.9) applied to (16.10) in place of (4.3): 

(16.18) \Ej\ < (l + Halloo + \\dp\\oo/co)Ei + 2^/Wi\\F[i\\ 
where cq is the constant in (1.6). By (16.14)-(16.17) we have 

r —1 

(16.19) 1177/11 < C (m(^_,||lT||, +p^_J|lT|U) +p^\\W\\r + ||F||, 

s=0 
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and using (16.13) 


r —1 


(16.20) 


\Hi\\ <K,Y^ {7hf_,+pf_,) (Cf + Ef) + K^p^{Cf + Ef) + ||F|| 


s=0 


Summing (16.18) over all I G 5 with |/| = r and using (16.20) we get 

. ^ I dEf 

(16.21) 


dt 


<Ki{l + ||(7||oo + ||5p||oo/co + Ell»Sp|U){C?+Ef) 




r — 1 


+ K,J2 i^f-s + Pr-s) {C^ + Ef) + l|i"l 


s=0 


Furthermore, by Lemma 12.1, (12.18) hold with U replaced by TZ and T replaced by S: 


(16.22) 


dC^ 


r— 1 


dt 


< 


{C^ + Ef) +K^Y. ^r-s (C'f + Ef) + ||F|| 


S=1 


(16.21) together with (16.22) gives us a bound for + E^ in terms of + Ef for s < r: 

(16.23) Cf (t) + Ef{t) < ’^^^(C^(O) + Ff (0)) 


i-t .’’“1 


+ K,e^^i:n<ir f (^Y^(^^n_^^pn_^^cf + Ef) + \\F\\r)dT 

s=l 


where n = 1 + H^Hoo + ||5p||oo/co + '^ses II^'S'pIIoo + ||<^||oo- Since we already have proven the bound 
for Eq = Fq in section 4, (16.23) inductively gives a bound for + E^. Hence by (16.13) we obtain: 


Lemma 16.1. Suppose that x,p G C'’'+^([0,r] x H), p\g^= 0, VwP 1^^^ < —cq < 0 anddivH = 0, where 
V = DfX. Let W be the solution of (16.1) where E is divergence free. Then there is a constant C 
depending only on the norm of {x,p), a lower hound for the constant cq and an upper bound for T, such 
that, for 0<t<T, we have 

(16.24) ||lF(t)||. + ||lH(t)||, + {W{t))A,r < C’(||1F(0)||, + ||VF(0)||, + (1F(0))a,. + ^ ||F||, dr) 
where 

(16.25) ||lH(t)||,= \\CijW{t)h.^n), {W{t))A,r= {E^sW{t), AC^sW{t))^/^ 

|7|<r, 7e7^ |/|<r.f65 


Note that the ||VF(t)||r is equivalent to the usual time independent Sobolev norm. Since there are 
compactly supported divergence free vector fields {W{t))A,r is only a semi-norm on divergence free 
vector helds, see (3.10). Furthermore, since 0 < cq < — Vjvp < C it follows from (3.11) that {W{t))A,r 
is equivalent to a time independent semi-norm given by (3.11) with / the distance function d{y), see 
(6.2). Since we only apply tangential vector fields, it also follows from (3.11) that, up to lower order 
terms that can be bounded by ||VF(t)||r., it is equivalent to that the normal component of the vector 
field Wn = NaW^ is in H^{dn). 

Definition 16.1. With notation as in (16.25) define F’'(H) to be the completion of C°°{Ll) in the norm 
||VF(t)||.r and define N'^{Q.) to be the completion of the divergence free C°°{Q.) vector fields in the norm 
\\WU. = \\W{t)\\r + {W{t))A,r. 

Since the projection onto divergence free vector fields is continuous in the norm it follows that 
id’’ is also the completion of the divergence free C°° vector fields in the F’’ norm. 
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Theorem 16.2. Suppose that x,p G C'’’'’"^([0, T] x fi), p\q^ = 0, V/vpl^^ < —cq < 0 and divDtx = 0. 
Then if initial data and the inhomogeneous term in (2.29) are divergence free and satisfy 

(16.26) {Wo,Wi) G X F G L^[0,T], {n)) 

the linearized equations (2.29) have a solution 

(16.27) (IT, IT) G C{[0,T],N^{n) x H^{n)) 

Proof. The existence of a solution in (16.27) follows from section 15 if initial data and the inhomogeneous 
term are divergence free and (7°° and the inhomogeneous term is supported in t > 0. By approximating 
the initial data and the inhomogeneous term in (16.26) with C°^ divergence free vector fields and 
applying the estimate (16.24) to the differences we get a convergent sequence in (16.27) so the limit 
must also be in this space. □ 
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